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ABSTRACT 

We show that as many as 4000 SNela may be required to detect the effect of weak 
lensing on their flux distribution with a high level of significance. However, if the 
intrinsic SNela magnitude dispersion is unknown one needs an even higher number of 
SNela (an order of magnitude more) to reach a similar level of statistical significance. 
Moreover, the ability to separate the lensing contribution from the intrinsic scatter 
depends sensitively on the amplitude of the latter. Using a Kolmogorov - Smirnov (K- 
S) test we check how the required number of SNela changes with level of significance. 
Our model incorporates a completely analytical description of weak lensing which has 
been tested extensively against numerical simulations. Thus, future missions such as 
SNAP may be able to detect non-Gaussianity at a lower significance level of 10% 
(through the K-S test) only if the intrinsic scatter is known from external data (e.g. 
from low redshift observations) whereas ALPACA with 100, 000 SNe will definitely 
detect non-Gaussianity with a very high confidence even if the intrinsic magnitude 
dispersion is not known a priori. 
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1 INTRODUCTION 

Type la supernovae (SNela) are powerful probes of the re- 
cent expansion of the universe and they provided the main 
contribution to the discovery of the present acceleration of 
the universe (Riess et al. 1998; Perlmutter et al. 1999). In- 
deed, SNela are standard candles with a small luminos- 
ity dispersion so that by measuring the flux received on 
the earth one can derive the luminosity distance of the 
source. Then, by observing many SNela one can measure 
the redshift-distance relation which provides constraints on 
cosmological parameters (Goobar & Perlmutter 1995). How- 
ever, even SNela are not perfect candles and are affected by 
various sources of noise such as the magnification produced 
by gravitational lensing, related to the fluctuations of the 
matter distribution along the line of sight (Kantowski et al. 
1995; Frieman 1997). Flux conservation implies that the ran- 
dom magnification shift is zero (Weinberg 1976) but weak 
lensing distortions increase the observed SNela magnitude 
scatter and lead to an extended high-luminosity tail (Wamb- 
sganss et al. 1997; Valageas 2000) . For a flux-limited survey 
weak lensing also leads to a slight bias towards larger lu- 
minosities close to the threshold (Valageas 2000) but this 
plays no significant role. Then, from the deviation of the 



magnitude distribution of 63 high redshift SNela (Riess et 
al. 2004) from a Gaussian, Wang (2005) claimed that weak 
lensing effects may have been detected. Although the distor- 
tion agrees at a qualitative level with the expectation from 
weak lensing magnification (i.e. there are three very bright 
SNela) the statistics was too small to draw firm conclusions. 
In this Letter we revisit this issue by investigating how many 
SNela are needed to detect with high confidence weak lens- 
ing effects. In sect. |2] we describe how weak gravitational 
lensing by large scale structures affects the apparent magni- 
tude of SNela. Next, assuming that the intrinsic magnitude 
fluctuation (including all sources of noise except lensing) is 
Gaussian with a known variance we discuss a Kolmogorov- 
Smirnov test to assess whether a sample of N supernovae 
may be drawn from such a Gaussian. Then, in sect. 13. ll wc 
compute for a survey such as the proposed SNAP 1 experi- 
ment (Aldering et al. 2004) at which A* a deviation from 
this Gaussian is detected with a high confidence level. In 
sect. 13.21 we generalize this procedure to the case where the 
intrinsic magnitude variance is unknown. We discuss the de- 
pendence of our results on the amplitude of the intrinsic 

1 http://snap.lbl.gov 
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SNela magnitude dispersion in sect. 13.31 and we conclude in 
sect.|U 



2 DETECTING WEAK LENSING 

If there are no distortions the apparent magnitude m app of a 
supernova at redshift z is related to its absolute magnitude 
Mt, and luminosity L„ by: 
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where <Il[z) is the luminosity distance, K(z) the "K- 
correction" which describes the redshift of the flux spectrum 
with respect to the observing filter and Lo the zero-point of 
the magnitude system. Therefore, inverting eq.Q observers 
who measure the flux from distant supernovae can derive 
cLl{z) which provides constraints on cosmological parame- 
ters. However, in practice one needs to take into account the 
intrinsic dispersion <7j nt of supernovae magnitudes, which is 
due to the dispersion of SNela luminosities themselves as 
well as to measurement noises and absorption along the line 
of sight. This implies a statistical analysis to extract the 
mean distance modulus m apP — A/*. Another distortion is 
due to weak lensing which can magnify the luminosity of dis- 
tant supernovae. Therefore, the apparent magnitude shows 
a fluctuation 5m around its average (m ap p}: 



5m — <5mi n t + 5m\ e 



(5m) = 0, 



(2) 



where we separate the intrinsic fluctuation Sm-mt and the 
gravitational lensing distortion ^ffMens- Note that since the 
mean weak lensing magnification is unity (as weak lensing 
only modifies the trajectory of light rays and does not change 
their energy) gravitational lensing does not bias the average 
apparent magnitude which allows one to derive dj, and to 
put constraints on cosmology. In this work, we focus on the 
fluctuating part 5m and we investigate how many SNela are 
required to detect weak lensing through the statistics of the 
fluctuations 5m, which depend on 5m\ cns . In the weak lens- 
ing regime which is appropriate for lensing by large scale 
structures that we consider here the magnification fj, is re- 
lated to the usual weak lensing convergence « by: 
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which can be written in terms of the density contrast 5 along 
the line of sight as: 
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Here Ho is the Hubble constant, \ ls the radial distance 
along the line of sight and T> the angular diameter distance. 
From the definition of magnitudes in Q we obtain for the 
apparent magnitude fluctuation: 



5m = 5mi n t + 5mi ell s = Smin 
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We shall assume in the following that 5mi n t is Gaussian 
with variance (Tint- Then, as in Valageas et al. (2005) the 
probability distribution of 5m can be written in terms of its 
generating function ips m as: 



P(5m) 
with: 
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Here we introduced the generating function 93<5m lens of the 
lensing fluctuation and we defined the ratio p of intrinsic 
and lensing variances by: 



P = 
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Finally, we use the model described in Barber et al. (2004) or 
Munshi et al. (2004) to obtain the generating functions of the 
convergence k whence of the lensing magnitude fluctuation 
<5miens, taking into account the redshift distribution of the 
sources (here SNela) as in Valageas et al. (2005). 

Then, from the observed distribution of apparent mag- 
nitudes, whence of 5m, one can recover the statistics of k. In 
this fashion, from the tails of the observed magnitude distri- 
bution Wang (2005) claimed that weak lensing effects may 
have been detected. However, the sample was too small (67 
high redshift SNela) to draw definite conclusions. Here we 
reconsider this question by computing how many SNela are 
needed to get a clear detection of weak lensing from SNela. 

To this order, we use a Kolmogorov-Smirnov (K-S) test 
(Press et al. 1986, Kendall & Stuart 1969) as follows. From 
a sample of N supernovae the observer can compare their 
magnitude distribution with a Gaussian Pg of variance ctg 
through the K-S distance d defined by: 



d — max \SN{5m) — Pg(< 5m)\ . 

dm 



(9) 



Here Pg(< 5m) is the trial cumulative Gaussian whereas 
SN(5m) is the discrete cumulative distribution obtained 
from the data. Thus, SN(Sm) is merely the fraction of ob- 
served SNela with a magnitude fluctuation smaller than 5m. 
As is well known, the interest of the K-S estimator d is 
that its distribution is universal when the data is compared 
with its parent distribution (null hypothesis) , whatever it is. 
Thus, the cumulative distribution of d writes in this case: 



P(> d) = Q KS {VNd) with QksW = -2^( 



(10) 



As expected, eg. (1101 shows that the distance d between the 
discrete cumulative distribution Sn and its continuous par- 
ent distribution scales as 1/yN. Then, by computing the 
distance d of the sample with respect to a trial Gaussian 
from eq.@, one can obtain from eg. (1101 the probability 
P(> d) that a distance of this size or larger would be ob- 
served if the trial Gaussian is the true parent magnitude 
distribution. Therefore, if this significance level P(> d) is 
smaller than a threshold P_ <<iC 1 one can conclude with good 
confidence that this trial Gaussian is not the true parent dis- 
tribution (disproof of the null hypothesis). In our case, this 
means that weak lensing has been detected since we assume 
that this is the only source of distortion from the Gaussian 
of intrinsic variance og with a a = a ln t- 

Thus, to find out for which N such a disproof of the 
Gaussian can be obtained with a high significance level we 
first choose a threshold P- <C 1 (for instance P- = 5%) 
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Figure 1. The probability distribution P(&m) (solid line) of the 
SNela magnitude fluctuation from the mean, as given by co.llrTl. 
The dashed line shows the Gaussian Pq of variance oq = (7i nt = 
0.1 mag which corresponds to neglecting weak- lensing effects. 

and compute from eq.QUJ the scaled variable A- such that 
Qks(A_) = P_ (for P_ = 5% this yields A_ = 1.34). Then, 
we draw a large number iV s i m of samples of N supernovae 
magnitude fluctuations Srrn (i = 1, .., N), for some value of 
N. As explained above, the distribution of these magnitudes 
is obtained from eq.©. Next, we compute for each sample k 
(k — 1, .., iVsi m ) the distance dk to the Gaussian of variance 
a a = ai n t, using eq.@. From this set {dk} we obtain the 
probability P(> A_) to measure a distance d larger than our 
threshold d_ = X-/s/N. This is simply the fraction of real- 
izations among our A?" s i m simulations with dk > d_. Then, we 
can repeat the same procedure for various N which provides 
the curve P(> A_ ; N) as a function of N (at fixed A_). Since 
the parent distribution © is different from the Gaussian of 
variance a- m t this probability P(> A_; N) increases with N 
and goes to unity at large N: for sufficiently large N we are 
sure to detect the difference between both PDF. Finally, we 
select a second threshold P+ ~ 1 (for instance P+ = 95%) 
and we find above which N, the probability P(> X-; N) 
becomes larger than P+. This value of TV* is the number of 
supernovae needed to detect with a high probability (P+) 
a weak lensing signature (defined as a distance d from the 
Gaussian which is more rare than P-). 



3 NUMERICAL RESULTS 

We assume a concordance ACDM cosmology with Q m = 
0.3, fi A = 0.7, a 8 = 0.88 and H = 70km/s/Mpc. We also 
adopt a redshift distribution of SNela as expected for the 
SNAP mission (Table 1 of Kim et al. 2004) which plans to 
observe about 6000 supernovae, of which 2000 may be used 
for cosmological purposes between redshifts of 0.1 and 1.7 
(Aldering et al. 2004). We also use throughout an intrinsic 
magnitude dispersion <7i nt =0.1 mag. 

We show in Fig.Qthe probability distribution P(Sm) of 
the SNela magnitude fluctuation 5m (solid line) from cq.Q. 
We also plot for comparison the Gaussian Pg of variance 
(Tint = 0.1. Thus, we see that weak-lensing effects increase 
the dispersion ((8m) 2 ) and distort the shape of the distri- 
bution with an extended bright tail. 
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Figure 2. The cumulative probability distribution P(> A) to 
measure a distance larger than d = A/ y/~N from the Gaussian of 
variance uq = cr; nt = 0.1 mag. The dot-dashed curves correspond 
to N = 100, 250, 500, 1000, 2000, 3000, 4000, 5000 from left 
to right. The left solid curve shows for reference the cumulative 
probability distribution of the distance from the parent distribu- 
tion JHJ. It is equal to Qks m eg. 1101 . 

3.1 Known intrinsic variance 

We first apply in this section the K-S test as described above 
in sect. [5] Thus, we show in Fig. [5] the cumulative prob- 
ability distribution P(> A) to measure a distance larger 
than d = \/y/N from the Gaussian of variance aa = <7i n t. 
These PDF are obtained for each N from the distribution 
of distances {dk} (k = 1, .., iV s i m ) associated with our N B im 
realizations of N supernovae. As iV increases the cumula- 
tive probability P(> A) develops a plateau which extends to 
larger values of A as it is easier to detect the deviation of the 
parent distribution © from the trial Gaussian of variance 
(Tint. Therefore, the probability P(> A_) grows with N. 

Thus, we show in Fig. |3 the curves P(> A_;7V) as a 
function of N, for three different thresholds P_ = 10%, 5% 
and 2.5% (corresponding to A_ = 1.20, 1.34 and 1.46) from 
top downto bottom. We can check that for low N the statis- 
tics is too small to obtain a clear detection of weak lensing 
and as iV increases the probability to measure the devia- 
tion from the Gaussian due to weak lensing effects grows 
to reach unity at N — > oo. Thus, we find that for signif- 
icance levels {P_ = 10%, P+ = 90%}, iV* = 2000 super- 
novae are sufficient to detect weak lensing. Higher levels 
{P_ = 5%,P+ = 95%} and {P_ = 2.5%, P+ = 97.5%} 
require N* = 3000 and N* = 4000 supernovae. 

3.2 Marginalizing over observed variance 

The procedure used in the previous section assumes that 
the intrinsic variance <Tj nt is exactly known so that any de- 
viation from the Gaussian of variance <Ti Ilt is interpreted as 
a detection of weak lensing. However, in practice the vari- 
ance (Ji n t is only known up to some finite accuracy. More- 
over, high redshift SNela may exhibit a somewhat different 
variance (because of the evolution of SNela metallicities, 
absorption by dust along the line of sight, etc.). Therefore, 
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Figure 3. The curves P(> A_ ; N) for three different thresholds 
P- = 10%, 5% and 2.5% (i.e. A_ = 1.20,1.34 and 1.46) from 
top downto bottom, as a function of N. At large N one detects 
almost surely (P(> A_;7V) ~ 1) a large deviation from the trial 
Gaussian (so large that it would have occurred with probability 
P— if the latter Gaussian had been the true parent distribution). 
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Figure 5. The curves P(> A_;7V) for three different thresholds 
P— = 10%, 5% and 2.5% from top downto bottom as in Fig. [3] 
but using the distance to the closest Gaussian displayed in Fig. HI 
A range of Gaussian PDFs were used (see text for more details). 
Triangles are actual estimates from our simulations whereas the 
solid lines are fitting functions. 
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Figure 4. Same as Fig. [5] but with marginalization over the 
observed variance. The dot-dashed curves show for various N 
the cumulative probability distribution P(> A m ; n ) to measure 
a distance larger than c( m i n = A m i n /-\/JV from the closest Gaus- 
sian among Gaussians of any variance (it is sufficient to span 
the range 0.09 < (tq < 0.14). They correspond to N = 
10000,20000,30000,40000,50000,60000 from left to right. The 
left solid curve shows for reference the distance from the parent 
distribution JHJ and obeys eg. 1101 . For each of these studies statis- 
tics are constructed from 1000 simulations. The intrinsic variance 
is again cr;, It = 0.1 mag. 

in this section we marginalize over the variance of the ob- 
served sample (keeping <7i nt = 0.1 mag for the true parent 
distribution), which implies that detection of weak lensing 
only depends on non-Gaussianities. Thus, for each realiza- 
tion k of N supernovae we compute all distances dk- p of this 
data set from an ensemble of trial Gaussians Pg. p of differ- 
ent variances (Jg-. p - From these dk- p we obtain the minimum 



distance d mln .k = min p {dk- p } . Thus d mln -k is the minimum 
distance between this realization and any Gaussian. In prac- 
tice we use a grid of variances <jg, p which spans the range 
[0.09, 0.14] with a step of 0.002. Obviously the distance dk- p 
increases at very small or very large variance aa and it is 
minimum for a G ~ of nt + of cns . For the cosmology and the 
redshift distribution that we use in this work we find that 
the minimum distance corresponds to ctg — 0.125. Then, 
from the distribution of minimum distances d m i n; ^ provided 
by our N 3 i m realizations of iV supernovae magnitudes we 
obtain the cumulative probability distribution P(> d m in) to 
observe a distance larger than d m i n from the closest possible 
Gaussian distribution. We show in Fig. [I] this cumulative 
probability distribution. Of course, we can check that for a 
given JV the typical distance d m i n is smaller than the dis- 
tance d to the fixed Gaussian of variance o G = Omt used 
in Fig. [5] Therefore, a larger number of SNela is needed to 
detect weak lensing. 

Applying the same procedure as in sect. l3.ll we can now 
display in Fig. E]the curves P(> A_;iV) as a function of N 
obtained from the minimum distance distributions shown 
in Fig. [1] We see that an order of magnitude more SNe are 
needed if the intrinsic variance is not known in advance (e.g. 
from low redshift studies). About 50,000 SNe are required 
to detect weak lensing effects in SNela studies through the 
K-S test with a high level of confidence for the redshift distri- 
bution that we have considered here. In particular, we now 
find that the significance levels {P_ = 10%, P+ = 90%}, 
{5%, 95%} and {2.5%, 97.5%} require JV» = 30,000,40,000 
and 45, 000 supernovae. Note that the SNAP mission actu- 
ally plans to observe ~ 10,000 SNela from which ~ 4000 
should be well-characterized (Albert et al. 2005). There- 
fore, it will be able to detect weak lensing only if the in- 
trinsic dispersion of SNela magnitudes is known. On the 
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Figure 6. The curves P(> A_ ; N) for three different thresholds 
P_ = 10%, 5% and 2.5% from top downto bottom. All simulations 
were performed with 3000 SNe. A set of 1000 simulations were 
performed to reduce the scatter. Solid lines represent a linear fit 
in the range of intrinsic variances considered in this study cr int = 
.1, .12, .14, .16. The variance of the trial Gaussian distribution is 
taken equal to the intrinsic variance cfq = <7j nt for this study 
(case of known intrinsic variance) . 

other hand, the JEDI 2 experiment proposes to observe over 
14, 000 SNela with well sampled light curve and good qual- 
ity spectra (Crotts et al. 2005) over < z < 1.7 whereas 
the ALPACA 3 experiment plans to observe > 100, 000 su- 
pernovae in the range 0.2 < z < 1 (Corasaniti et al. 2005). 

3.3 Dependence on intrinsic variance 

In these studies we have assumed that the intrinsic variance 
(which can be unknown) is (Tint = 0.1 mag. However these re- 
sults are quite sensitive to (Tint- In Fig.HJwe plot for the case 
of N = 3000 SNela the cumulative probability P(> A_; N) 
as a function of (Tint- We consider the three thresholds used 
in Figs.|3]|5] and we use 1000 simulations. As in sect. l3.ll we 
consider the case where the intrinsic variance is known so 
that the observed SNela sample is compared with the trial 
Gaussian of variance gg = (Tint- Of course, for low (Tint it 
is easy to detect weak lensing (the probability P(> A-; N) 
goes to unity) since the amplitude of weak lensing effects 
becomes larger than the intrinsic dispersion of SNela magni- 
tudes whereas for high <7j nt weak lensing distortions become 
relatively negligible (P(> X~;N) goes to zero). We can see 
that this probability P(> A_; N) is quite sensitive to (Tint as 
it exhibits a fast decrease for larger (Tint. This implies that 
the number of SNela required to detect weak lensing signa- 
tures through the Kolmogorov-Smirnov test grows quickly 
with the intrinsic SNela magnitude variance. In particular, 
for crint = 0.16 mag we find that 7000, 10000 and 20000 
SNela are required in order to achieve the confidence levels 
{10%, 90%}, {5%, 95%} and {2.5%, 97.5%} (in the case of 
known intrinsic variance as in sect. I3.lt . 

2 http://jedi.nhn.ou.edu 

3 http:/ /www. astro. ubc.ca/LMT/alpaca/index. html 



4 CONCLUSIONS AND OUTLOOK 

In this Letter we have addressed the issue of determining 
the number of observed SNela beyond which weak lensing 
effects can be detected with a high confidence. For the con- 
cordance ACDM cosmology, using a model of the large-scale 
matter distribution which has been checked against numer- 
ical simulations, we found that 4000 SNela are necessary 
to distinguish a weak lensing signature with a significance 
level of 2.5% through a Kolmogorov-Smirnov test. To reach a 
significance level of 10% we only need 2000 SNela. This pro- 
cedure compares the magnitude distribution of the observed 
SNela with a Gaussian of fixed variance, assuming that the 
latter describes all sources of noise except for weak lensing 
magnification. If we consider the variance to be a free pa- 
rameter (e.g. the intrinsic SNela magnitude dispersion or the 
instrumental noise are not accurately known beforehand) we 
find that 45, 000 supernovae are required to detect with a 
high confidence (at a 2.5% level) non-Gaussian signatures. 
Therefore, future experiments such as those planned within 
the Joint Dark Energy Mission will exhibit clear weak lens- 
ing signatures if the intrinsic magnitude dispersion of SNela 
is well known. However to be more confident without any a 
priori knowledge of (Tint we will have to wait for surveys such 
as ALPACA. Of course, the possibility of detecting weak 
lensing effects on SNela magnitude distributions also im- 
plies that such gravitational lensing effects should be taken 
into account or used as a complementary tool to constrain 
cosmology (e.g., Dodelson & Vallinotto 2005). 
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